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A new type of self-similar potential is used to study a multibarrier system made of graphene. 
Such potential is based on the traditional middle third Cantor set rule combined with a scaling of 
the barriers height. The resulting transmission coefficient for charge carriers, obtained using the 
quantum relativistic Dirac equation, shows a surprising self-similar structure. The same potential 
does not lead to a self-similar transmission when applied to the typical semiconductors described 
by the non-relativistic Schrodinger equation. The proposed system is one of the few examples in 
which a self-similar structure produces the same pattern in a physical property. The resulting 
scaling properties are investigated as a function of three parameters: the height of the main barrier, 
the total length of the system and the generation number of the potential. These scaling properties 
are first identified individually and then combined to find general analytic scaling expressions. 
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GraphenejlJ is considered as one of the most promising new materials. This first truly two-dimensional crystal 0 
has impressive physical properties 00 - As a consequence, graphene is not only important from a technological point 
of view, but also is considered as an inflexion point in quantum physics. For instance, charge carriers in graphene 
follow an effective quantum relativistic (Dirac) equation instead of the usual Schrodinger one. This leads to new 
effects 0-0. Specially, the scattering produced by external potentials or impurities is different from what is observed 
in ordinary materials. Charge carriers can be transmitted perfectly due to the Klein effect and random sustitutional 
impurities can produce multifractal states 00. Thus, a whole new world is open to investigate the possible effects 
of external potential geometries, which can be imposed to graphene by different means, as for example, using certain 
type of substrates, strain, electrostatic gates, impurities, electromagnetic fields, etc. [1. Ho|. 


For applications, an important issue is how to engineer substrate-induced bandgap opening in epitaxial graphene 
M- This leads to investigate how different multibarrier geometries affect the properties of graphene. Among these 
possible geometries, self-similarity has a paramount importance, since scale invariance is a fundamental property of 
many natural phenomena, as can be corroborated in reports that ranges from the distribution and abundance of 
species, temporal occurrence of earthquakes, and even in the growth of complex networks and trees mm. 

From the technological standpoint, self-similarity can also be exploited to produce useful devices wEa . Par¬ 
ticularly, photoconductive fractal antennas show an efficient multiband emission of terahertz radiation owing to the 
self-similarity of the fractal structure [16!|. 

Despite that self-similar structures are getting plenty of attention (Ts) , only a handful of experiments corroborate a 
self-similar behavior in its physical properties [19|-[21]. Likewise, many theoretical works claim that physical properties, 
such as tranmission or reflection probabilities, display self-similarity. Unfortunately, most of them are just a matter 
of visual perception, since the scaling properties are never reported. Even in the case of quasicrystals, which are 
considered as one of the archetypical examples of self-similar structures, it has been elusive to find the scaling 
properties of the corresponding physical properties 02|, |23[ . 

Here we propose a novel self-similar multibarrier structure in graphene in which the main physical property, namely 
charge transport, presents a surprising self-similar pattern. This self-similar structure does not produce a similar 
behavior in a usual semiconducting material described by the Schrodinger equation. Thus, two essential ingredients 
are shown to be into play: a self-similar structure and quantum relativistic mechanics. 

To build our multibarrrier structure in one direction (the x axis), we implement a variant of the middle third Cantor 
case [24j (see supplementary material) using square barriers as indicated in Fig. [TJ We begin with a line-segment of 
length Lt, we divide it in thirds and put a rectangular barrier of energy Vo in the middle third (its width is Lt/ 3). 
We call this the generation one (Gl). The second step consists in taking the remaining line-segments, divide them 
into thirds and place within scaled copies of the existing barrier. We scale the height and width of the barrier by 
a scaling factor of 1/3. Additionally the existing barrier is also scaled in its width. Since scaling is a contraction 
transformation, each scaled copy will have a unique fixed point, see the fixed point theorem in metric spaces [25j . 
By fixed point we meant an element of the function’s domain that is mapped to itself by the function. In this case, 
this point is a fixed side of the barriers, which will be chosen as the right-hand side. So, every successive iteration 
(generation) ( Gj ) of the structure is obtained by repeating the second step of the construction for the remaining 
line-segmentes, using as the “existing barrier” those added in the last iteration, as shown in Fig. |1] 


This potential can be described as follows. The j- th generation of the potential will consists of N G j barriers; the 
corresponding barriers (those added in the j-th iteration) has a width W G (j ) and energy (height) V G j. Such barrier 
is obtained as a scaled replica of barriers from generation G(j — 1), 

N Gj = 2 Gi - 1, 

w Lt 

W G j = a = -TT- 


Vg,' = 


3 GU 

^0 

3 Gj-i' 


This bidimensional potential does not depend on the perpendicular direction (y axis), i.e., V(x,y) = V{x) and can 
be adapted to graphene 0 2fil-l28l| by means of symmetry breaking substrates [2fil - l3i| . by gated graphene or in other 
systems, like optical analogies, with a linear spectrum |32 - l34l |. 

Let us now study the transmission coefficient of charge carriers in this multibarrier system. For graphene, we use 
the Viana-Gomes et al. formalism [30j to compute the transmission coefficient of electrons in gapped graphene. The 
Hamiltonian for the Dirac electrons under a potential V ( x ) is given by 


H — v F a -p + a z V(x), (1) 

where V(x) = m(x)vp ii x £ [xi,X 2 ] and zero in other case. V(x) defines the region where there is a mass term 
(q- region) and where carriers are massles Dirac fermions (fc-region). The first term of the Hamiltonian is just the 
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FIG. 1: Potential construction. First four generations (G1,G2,G3,G4) of the multibarrier potential. Each 
barrier is scaled (we denote the scaling factor by a) in its width, and so the barrier width is Lt/a. The barriers 
added in the last iteration are scaled in its height. For example, the height of the green ones are one third of the 
height of the blue ones. The zooms in G3 and G4 illustrate how the potential resembles a previous generation after 

an appropiate rescaling. 


Dirac equation used to describe electrons in graphene, where vf is the Fermi velocity, a the set of Pauli matrices and 
p the moment. For normal incidence on the x direction the corresponding wavefunctions are [3fj 

= 71 (-1 )'“•'■ < 2 > 

and 

^ W = 75 ( „1) < 3 > 

where u± = ±sign(U) and v± = hV p(^ q ^\ k ) • Using these equations, one establishes the continuity conditions when 
the wave goes from the fc-region to the g-region. 


With the potential described in Fig. [Qat hand, we compute the transmission coefficient, defined as 


T G (E) 


l^(Q)l 2 

\^k( Lt W 


( 4 ) 


using the transfer matrix formalism [30j, l35i . [36 ]. The transfer matrix can be obtained by imposing the continuity 
conditions of the wave functions in each region of the self-similar structure. From this we can extract the transmission 
coefficient for a given energy, see for instance [30j, [H]. 

Since the construction of the multibarrier structure is based on the length-scaling of the previous generation barriers, 
and the addition of energy-scaled new barriers, we search for relations between transmission curves corresponding to 
different generations of the multibarrier structure. So, we fix the total length of the structure, denoted by Lt (in this 
case we use the value of Lt = 3500 A) and the height of the main barrier, denoted by Vq (using the value of Vq = 1 eV). 
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c) Comparison, Schrodinger case. d) Scaling, Schrodinger case. 



FIG. 2: Transmission coefficient scaling for different generations, a) Transmission as a function of energy 
for the sixth (G = 6) and seventh (G = 7) generation of multibarrier graphene using the parameters Vo = 1.0 eV and 
Lt = 3500 A. b) Scaled transmission for G = 6 and G = 7. The scaling transformation is given by the power law Eq. 
ED, which gives an excellent agreement of the curves. In c) and d), the transmission and scaled transmission 
respectively are presented for the same potential but for a typical semiconductor which follows the Schrodinger 

equation. No self-similar scaling is observed. 


In figure[2Di we show the transmission curves corresponding to generations 6 and 7 of the potential, denoted by Tq(E) 
and Tf(E) respectively. As we can see, it is qualitatively clear that they are related by some kind of transformation. 
This transformation turns out to be Tq(E) = [Tt(E )] 9 , where the subindex stands for the generation number, as seen 
in hgure[5Jr. Thus, T 6 (E) are Ty(E) scaled by a single exponent. As may be expected, other bigger generations also 
follow the same rule. More formaly 

T g (E) « [T G+1 (E)} 9 , (5) 

where G > 6 and G stands for the generation number. This means that any two transmission curves that differ 
only in the generation of their correspondig potential, and are related by some power of the form 9 m , where m is 
the difference between the generations. This suggests that the transmission curves follows an scaling rule which is a 
signature of the self-similar phenomena. In fact, here we are in presence of the first true self-similar physical property 
of graphene. Morever, even for perfect fractals it is known that the physical properties usually follow a multifractal 
behavior [23|, [37] but here only one single fractal exponent is obtained. 

It is important to remark here that there are always finite size effects due to the break in the scaling of the potential, 
which affects the first generations. In this particular case, one cannot extend self-similarity beyond the length of the 
system. This effect is explained in the supplementary section where we present the transmittance from generation 
G1 to G9. For low generations, the transmittance is basically dominated by the biggest energy barrier. Between G5 
and G6, there is a dramatic change in the behavior since now the system has transmittance for energies lower than 
the main barrier height. As generations increase, scaling appears. Generation 6 is the crossover since as shown in 
the supplementary section, geometrical finite size effects in the potential die out after it. Also, is clear that in real 
systems there is also another finite size effect at the bottom of the hierarchy, where the Dirac approach can not be 
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made due to the rapid variation of the potential at atomic distances. 

To highlight the difference with a typical Schrodinger semiconductor, in figure [21 panels c) and d), we present the 
transmission calculated for the same potential generations using Al^Gai-^As/GaAs. In such system, the height of 
the quantum wells are controlled by the Aluminium concentration ( x ). 

So far, we have found a scaling between transmission curves that correspond to different generations of the potential, 
but it seems possible to extend the search for scaling features to other parameters. We have two options at hand, the 
energy of the starting barrier and the total length of the structure. 


a) Changing the main barrier energy. b) Scaling in the main barrier energy. 



c) Changing the total length. d) Scaling in the total length. 



FIG. 3: Scaling between barrier energy, a) Comparison between transmission curves for different values of Vo 
using generation 6. This case corresponds to the transmission for two potentials with different heights of the main 
barrier, in this case Vo = 0.5 eV and Vo = 1.0 eV. All other parameters are equal, b) Comparison between the scaled 
plots of a). The scaling consists in raising Vo = 0.5 eV to the fourth power. Scaling between lengths, c) 
Comparison between the transmission curves for different lengths (Lt). In this case it is quite clear that the scaling 
cannot be just raising one curve to a certain power, d) Scaling of the curve with length Lt = 500 A. This scaling is 
made of two parts, raising the curve to the fourth power and multiplying its energy by two. 


Let us explore now the first option, that is, we take a fixed generation of the potential, say generation 6, and make 
the calculation for two starting energies. Lets take Vo = 0.5 eV and Vo = 1.0 eV. In this case, the appearance of the 
curves seen in figure [3jr, resembles those corresponding to a change in the generation of the potential. In the case 
of the nonrelativistic transmission coefficient (e. g. GaAs), for the same geometric set up, one observes in general 
a shift of the first transmission resonance. Nevertheless, in this case the general behavior of the transmission curve 
is qualitatively the same. Due to the resemblance of these curves with the ones of changing generations, we try the 
same kind of transformation but with a different power. In this case, one has to raise the curve corresponding to 
Vo = 0.5 eV to the fourth power to get the one corresponding to Vo = 1.0 eV, see figure[3jD. Thus T\,q(E) r; T 0 . 5 (E) a , 
where the subindex stands for the Vo of the respective structure. In fact, this scaling behavior is more general, i. e., 
it applies for every multiple of two, resulting in the expression 

T iVo (Ef *T Vo (E). (6) 

Finally we compare curves corresponding to different lengths of the whole structure, denoted by Lt. Although 
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nothing suggests scaling in this variant at all, see figure ©:, our previous results lead us to a different conclusion. 
From the results for the case of changing Vo, we try to find a similar scaling law for Lt, but in this case, the scaling 
in the energy gives the relation, 

T XLt (-E) a2 ^T Lt (E). ( 7 ) 

“ a 

With this scaling rule applied to the curve corresponding to 500 A, a very good approximation is obtained, see figure 

© 1 . 

Combining expressions ©, © and ED, it is possible to get a good approximation to almost every curve that results 
from the combination of parameters. To do this, all the parameters and the variable E will be arguments of T, 

T(E, G, V 0 , Lt) « T(aE, G - m,\v 0 , aLt) k2/a29m , (8) 

k 

for G — m > 6. All these approximations are quite useful to predict the transmission as function of the parameters, 
but it is necesary to test the correctness. So we calculate the root mean square of the difference between the target 
curve and the scaled one, we call it drms. For parameter G (scaling between generations of the potential, see figure 
|2b) the drms is 5.50 x 10~ 1 2 3 4 5 6 7 ; for Vq (scaling in the energy of the main barrier, see figure©)) the drms is 2.05 x 10 -4 ; 
at last for Lt (scaling in the total length of the multibarrier structure, see figure©) the drms is 1.43 x 10 -3 . From 
this it is quite clear that the better approximation is for the scaling between generations. 

This lead us to examinate a wider range of energy and check our scaling-laws values (see supplementary material). 
An analysis of the results, shows that the fit is excellent for small energies but tends to decrease in quality as the 
energy increases, as one should expect in the limiting case of very high energies. Also, the error tends to grow for 
small transmittance energies. However, as seen in the supplementary material, if the error is computed for bigger 
generations, for example, for G = 7 against scaled G = 8, then the error decreases in a dramatic way, as expected for 
a good scaling relation. The drms for the Schroedinger case does not reduce as the generations grow, suggesting the 
absence of a fixed point. 

The numerical results presented above suggests that this type of potential indeed shows scaling properties, although 
this scaling is not found in a single curve, but between curves corresponding to different parameters. This is a very 
interesting result because in fact the transmission coefficient is an intricate combination of the wave functions, and 
even more, for a multibarrier potential it is the result of the product of several transfer matrices. This s ugg ests that 
the wave function have self-similar properties. That is the affirmation in the Gaggero-Pujals theorem [381 ] for non- 
relativistic electrons. Thus, it seems plausible that an extension of that theorem could be established for relativistic 
equations. Another important possibility is to study the nonperpedicular incidence, since Klein tunneling is ruled out 
in this case [39j . 

In conclusion, the self-similar potential proposed here show various kinds of scaling properties in the transmission 
curves. This scaling properties gave the possibility to obtain, given a transmission curve and its set of parameters, an 
almost perfect approximation of the curve correspondig to certain transformation of the original parameters, see ©. 
It is also suggested that the scaling properties of the transmission curves are drastically affected by the mix of two 
geometric symmetries, for example, reflection and self-similarity, see 0, [40| . We believe that the proposed potential 
is ideal to study scaling properties in multibarrier systems, although our preliminary results considering electrostatic 
self-similar barriers in grafene indicate also the presence of a self-similar behavior. 

Finally, we found that two ingredients are needed in order to obtain a pure monofractal scaling of the transmittance, 
a self-similar barrier and quantum relativistc equations , since the non-relativistic version does not display scaling. A 
possible heuristic explanation for such phenomena lies in the linear energy-momentum dispersion relation of the Dirac 
equation. 
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